Data-specific substitution models improve protein-based phylogenetics.

Calculating amino-acid substitution models that are specific for individual protein data sets is often difficult due to the computational burden of estimating large numbers of rate parameters. In this study, we tested the computational efficiency and accuracy of five methods used to estimate substitution models, namely Codeml, FastMG, IQ-TREE, P4 (maximum likelihood), and P4 (Bayesian inference). Data-specific substitution models were estimated from simulated alignments (with different lengths) that were generated from a known simulation model and simulation tree. Each of the resulting data-specific substitution models was used to calculate the maximum likelihood score of the simulation tree and simulated data that was used to calculate the model, and compared with the maximum likelihood scores of the known simulation model and simulation tree on the same simulated data. Additionally, the commonly-used empirical models, cpREV and WAG, were assessed similarly. Data-specific models performed better than the empirical models, which under-fitted the simulated alignments, had the highest difference to the simulation model maximum-likelihood score, clustered further from the simulation model in principal component analysis ordination, and inferred less accurate trees. Data-specific models and the simulation model shared statistically indistinguishable maximum-likelihood scores, indicating that the five methods were reasonably accurate at estimating substitution models by this measure. Nevertheless, tree statistics showed differences between optimal maximum likelihood trees. Unlike other model estimating methods, trees inferred using data-specific models generated with IQ-TREE and P4 (maximum likelihood) were not significantly different from the trees derived from the simulation model in each analysis, indicating that these two methods alone were the most accurate at estimating data-specific models. To show the benefits of using data-specific protein models several published data sets were reanalysed using IQ-TREE-estimated models. These newly estimated models were a better fit to the data than the empirical models that were used by the original authors, often inferred longer trees, and resulted in different tree topologies in more than half of the re-analysed data sets. The results of this study show that software availability and high computation burden are not limitations to generating better-fitting data-specific amino-acid substitution models for phylogenetic analyses.

Recoding Amino Acids to a Reduced Alphabet may Increase or Decrease Phylogenetic Accuracy.

Common molecular phylogenetic characteristics such as long branches and compositional heterogeneity can be problematic for phylogenetic reconstruction when using amino acid data. Recoding alignments to reduced alphabets before phylogenetic analysis has often been used both to explore and potentially decrease the effect of such problems. We tested the effectiveness of this strategy on topological accuracy using simulated data on four-taxon trees. We simulated alignments in phylogenetically challenging ways to test the phylogenetic accuracy of analyses using various recoding strategies together with commonly used homogeneous models. We tested three recoding methods based on amino acid exchangeability, and another recoding method based on lowering the compositional heterogeneity among alignment sequences as measured by the Chi-squared statistic. Our simulation results show that on trees with long branches where sequences approach saturation, accuracy was not greatly affected by exchangeability-based recodings, but Chi-squared-based recoding decreased accuracy. We then simulated sequences with different kinds of compositional heterogeneity over the tree. Recoding often increased accuracy on such alignments. Exchangeability-based recoding was rarely worse than not recoding, and often considerably better. Recoding based on lowering the Chi-squared value improved accuracy in some cases but not in others, suggesting that low compositional heterogeneity by itself is not sufficient to increase accuracy in the analysis of these alignments. We also simulated alignments using site-specific amino acid profiles, making sequences that had compositional heterogeneity over alignment sites. Exchangeability-based recoding coupled with site-homogeneous models had poor accuracy for these data sets but Chi-squared-based recoding on these alignments increased accuracy. We then simulated data sets that were compositionally both site- and tree-heterogeneous, like many real data sets. The effect on the accuracy of recoding such doubly problematic data sets varied widely, depending on the type of compositional tree heterogeneity and on the recoding scheme. Interestingly, analysis of unrecoded compositionally heterogeneous alignments with the NDCH or CAT models was generally more accurate than homogeneous analysis, whether recoded or not. Overall, our results suggest that making trees for recoded amino acid data sets can be useful, but they need to be interpreted cautiously as part of a more comprehensive analysis. The use of better-fitting models like NDCH and CAT, which directly account for the patterns in the data, may offer a more promising long-term solution for analyzing empirical data. [Compositional heterogeneity; models of evolution; phylogenetic methods; recoding amino acid data sets.].

Integrated phylogenomic analyses unveil reticulate evolution in Parthenocissus (Vitaceae), highlighting speciation dynamics in the Himalayan-Hengduan Mountains.

Hybridization caused by frequent environmental changes can lead both to species diversification (speciation) and to speciation reversal (despeciation), but the latter has rarely been demonstrated. Parthenocissus, a genus with its trifoliolate lineage in the Himalayan-Hengduan Mountains (HHM) region showing perplexing phylogenetic relationships, provides an opportunity for investigating speciation dynamics based on integrated evidence. We investigated phylogenetic discordance and reticulate evolution in Parthenocissus based on rigorous analyses of plastome and transcriptome data. We focused on reticulations in the trifoliolate lineage in the HHM region using a population-level genome resequencing dataset, incorporating evidence from morphology, distribution, and elevation. Comprehensive analyses confirmed multiple introgressions within Parthenocissus in a robust temporal-spatial framework. Around the HHM region, at least three hybridization hot spots were identified, one of which showed evidence of ongoing speciation reversal. We present a solid case study using an integrative methodological approach to investigate reticulate evolutionary history and its underlying mechanisms in plants. It demonstrates an example of speciation reversal through frequent hybridizations in the HHM region, which provides new perspectives on speciation dynamics in mountainous areas with strong topographic and environmental heterogeneity.

Comparative mitogenomic analyses and gene rearrangements reject the alleged polyphyly of a bivalve genus.

Background: The order and orientation of genes encoded by animal mitogenomes are typically conserved, although there is increasing evidence of multiple rearrangements among mollusks. The mitogenome from a Brazilian brown mussel (hereafter named B1) classified as Perna perna Linnaeus, 1758 and assembled from Illumina short-length reads revealed an unusual gene order very different from other congeneric species. Previous mitogenomic analyses based on the Brazilian specimen and other Mytilidae suggested the polyphyly of the genus Perna. Methods: To confirm the proposed gene rearrangements, we sequenced a second Brazilian P. perna specimen using the "primer-walking" method and performed the assembly using as reference Perna canaliculus. This time-consuming sequencing method is highly effective when assessing gene order because it relies on sequentially-determined, overlapping fragments. We also sequenced the mitogenomes of eastern and southwestern South African P. perna lineages to analyze the existence of putative intraspecific gene order changes as the two lineages show overlapping distributions but do not exhibit a sister relationship. Results: The three P. perna mitogenomes sequenced in this study exhibit the same gene order as the reference. CREx, a software that heuristically determines rearrangement scenarios, identified numerous gene order changes between B1 and our P. perna mitogenomes, rejecting the previously proposed gene order for the species. Our results validate the monophyly of the genus Perna and indicate a misidentification of B1.

A low-latitude species pump: Peripheral isolation, parapatric speciation and mating-system evolution converge in a marine radiation.

Geologically recent radiations can shed light on speciation processes, but incomplete lineage sorting and introgressive gene flow render accurate evolutionary reconstruction and interpretation challenging. Independently evolving metapopulations of low dispersal taxa may provide an additional level of phylogeographic information, given sufficiently broad sampling and genome-wide sequencing. Evolution in the marine brown algal genus Fucus in the south-eastern North Atlantic was shaped by Quaternary climate-driven range shifts. Over this timescale, divergence and speciation occurred against a background of expansion-contraction cycles from multiple refugia, together with mating-system shifts from outcrossing (dioecy) to selfing hermaphroditism. We tested the hypothesis that peripheral isolation of range edge (dioecious) F. vesiculosus led to parapatric speciation and radiation of hermaphrodite lineages. Species tree methods using 876 single-copy nuclear genes and extensive geographic coverage produced conflicting topologies with respect to geographic clades of F. vesiculosus. All methods, however, revealed a new and early diverging hermaphrodite species, Fucus macroguiryi sp. nov. Both the multispecies coalescent and polymorphism-aware models (in contrast to concatenation) support sequential paraphyly in F. vesiculosus resulting from distinct evolutionary processes. Our results support (1) peripheral isolation of the southern F. vesiculosus clade prior to parapatric speciation and radiation of hermaphrodite lineages-a "low-latitude species pump". (2) Directional introgressive gene flow into F. vesiculosus around the present-day secondary contact zone (sympatric-allopatric boundary) between dioecious/hermaphrodite lineages as hermaphrodites expanded northwards, supported by concordance analysis and statistical tests of introgression. (3) Species boundaries in the extensive sympatric range are probably maintained by reproductive system (selfing in hermaphrodites) and reinforcement.

metaGOflow: a workflow for the analysis of marine Genomic Observatories shotgun metagenomics data.

BACKGROUND: Genomic Observatories (GOs) are sites of long-term scientific study that undertake regular assessments of the genomic biodiversity. The European Marine Omics Biodiversity Observation Network (EMO BON) is a network of GOs that conduct regular biological community samplings to generate environmental and metagenomic data of microbial communities from designated marine stations around Europe. The development of an effective workflow is essential for the analysis of the EMO BON metagenomic data in a timely and reproducible manner. FINDINGS: Based on the established MGnify resource, we developed metaGOflow. metaGOflow supports the fast inference of taxonomic profiles from GO-derived data based on ribosomal RNA genes and their functional annotation using the raw reads. Thanks to the Research Object Crate packaging, relevant metadata about the sample under study, and the details of the bioinformatics analysis it has been subjected to, are inherited to the data product while its modular implementation allows running the workflow partially. The analysis of 2 EMO BON samples and 1 Tara Oceans sample was performed as a use case. CONCLUSIONS: metaGOflow is an efficient and robust workflow that scales to the needs of projects producing big metagenomic data such as EMO BON. It highlights how containerization technologies along with modern workflow languages and metadata package approaches can support the needs of researchers when dealing with ever-increasing volumes of biological data. Despite being initially oriented to address the needs of EMO BON, metaGOflow is a flexible and easy-to-use workflow that can be broadly used for one-sample-at-a-time analysis of shotgun metagenomics data.

Biochemical and molecular responses of the Mediterranean mussel (Mytilus galloprovincialis) to short-term exposure to three commonly prescribed drugs.

Pharmaceuticals represent a group of emerging contaminants. The short-term effect (3 and 7 days) of warfarin (1 and 10 mg L-1), dexamethasone (0.392 and 3.92 mg L-1) and imidazole (0.013 and 0.13 mg L-1) exposure was evaluated on mussels (Mytilus galloprovincialis). Total antioxidant status, glutathione reductase, glutathione peroxidase (GPx) and superoxide dismutase enzyme activities, and the expression of genes involved in the xenobiotic response (ATP binding cassette subfamily B member 1 (abcb1) and several nuclear receptor family J (nr1j) isoforms), were evaluated. All nr1j isoforms are suggested to be the xenobiotic receptor orthologs of the NR1I family. All drugs increased GPx activity and altered the expression of particular nr1j isoforms. Dexamethasone exposure also decreased abcb1 expression. These findings raised some concerns regarding the release of these pharmaceuticals into the aquatic environment. Thus, further studies might be needed to perform an accurate environmental risk assessment of these 3 poorly studied drugs.

Inferring the Deep Past from Molecular Data.

There is an expectation that analyses of molecular sequences might be able to distinguish between alternative hypotheses for ancient relationships, but the phylogenetic methods used and types of data analyzed are of critical importance in any attempt to recover historical signal. Here, we discuss some common issues that can influence the topology of trees obtained when using overly simple models to analyze molecular data that often display complicated patterns of sequence heterogeneity. To illustrate our discussion, we have used three examples of inferred relationships which have changed radically as models and methods of analysis have improved. In two of these examples, the sister-group relationship between thermophilic Thermus and mesophilic Deinococcus, and the position of long-branch Microsporidia among eukaryotes, we show that recovering what is now generally considered to be the correct tree is critically dependent on the fit between model and data. In the third example, the position of eukaryotes in the tree of life, the hypothesis that is currently supported by the best available methods is fundamentally different from the classical view of relationships between major cellular domains. Since heterogeneity appears to be pervasive and varied among all molecular sequence data, and even the best available models can still struggle to deal with some problems, the issues we discuss are generally relevant to phylogenetic analyses. It remains essential to maintain a critical attitude to all trees as hypotheses of relationship that may change with more data and better methods.

Author Correction: Phylogenomics provides robust support for a two-domains tree of life.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Uncovering the shell game with barcodes: diversity of meiofaunal Caecidae snails (Truncatelloidea, Caenogastropoda) from Central America.

Caecidae is a species-rich family of microsnails with a worldwide distribution. Typical for many groups of gastropods, caecid taxonomy is largely based on overt shell characters. However, identification of species using shell characteristics is problematic due to their rather uniform, tubular shells, the presence of different growth stages, and a high degree of intraspecific variability. In the present study, a first integrative approach to caecid taxonomy is provided using light-microscopic investigation with microsculptural analyses and multi-marker barcoding, in conjunction with molecular species delineation analyses (ABGD, haplotype networks, GMYC, and bPTP). In total 132 specimens of Caecum and Meioceras collected during several sampling trips to Central America were analyzed and delineated into a minimum of 19 species to discuss putative synonyms, and supplement the original descriptions. Molecular phylogenetic analyses suggest Meioceras nitidum and M. cubitatum should be reclassified as Caecum, and the genus Meioceras might present a junior synonym of Caecum. Meiofaunal caecids morphologically resembling C. glabrum from the Northeast Atlantic are a complex of cryptic species with independent evolutionary origins, likely associated with multiple habitat shifts to the mesopsammic environment. Caecum invisibile Egger & Jörger, sp. nov. is formally described based on molecular diagnostic characters. This first integrative approach towards the taxonomy of Caecidae increases the known diversity, reveals the need for a reclassification of the genus Caecum and serves as a starting point for a barcoding library of the family, thereby enabling further reliable identifications of these taxonomically challenging microsnails in future studies.

Prokaryotic diversity in stream sediments affected by acid mine drainage.

The microbial communities in mining impacted areas rely on a variety of mechanisms to survive in such extreme environments. In this work, a meta-taxonomic approach using 16S rRNA gene sequences was used to investigate the prokaryotic diversity of sediment samples from water bodies affected by acid mine drainage at the São Domingos mining area in the south of Portugal. Samples were collected in summer and winter from the most contaminated sites from where the water flows downstream to the freshwater of Chança's river reservoir. The prokaryotic diversity on water bodies' sediments allowed us to distinguish the highly contaminated sites (pH ≈ 2) from sites with intermediate levels of contamination (pH ≈ 3-6.5), and from sites without contamination (pH ≈ 7.5). The abundances of acidophiles of genera Acidiphilium, Acidibacter, Acidobacterium and Acidocella in the sediments were correlated with the level of acid mine drainage contamination. The two first genera were among the 30 most abundant prokaryotes in all contaminated samples, including one (SS2w), where the contamination was very diluted, thereby emphasizing the impact that such type of pollution can have in the microbial communities of sediments. In addition, the high abundances of archaeal taxa from class Thermoplasmata and of bacteria from family RCP1-48 in the sediments from the most contaminated site corroborate their importance in such ecosystems and a putative role in the generation of acid mine drainage.

The Chloroplast Land Plant Phylogeny: Analyses Employing Better-Fitting Tree- and Site-Heterogeneous Composition Models.

The colonization of land by descendants of charophyte green algae marked a turning point in Earth history that enabled the development of the diverse terrestrial ecosystems we see today. Early land plants diversified into three gametophyte-dominant lineages, namely the hornworts, liverworts, and mosses, collectively known as bryophytes, and a sporophyte-dominant lineage, the vascular plants, or tracheophytes. In recent decades, the prevailing view of evolutionary relationships among these four lineages has been that the tracheophytes were derived from a bryophyte ancestor. However, recent phylogenetic evidence has suggested that bryophytes are monophyletic, and thus that the first split among land plants gave rise to the lineages that today we recognize as the bryophytes and tracheophytes. We present a phylogenetic analysis of chloroplast protein-coding data that also supports the monophyly of bryophytes. This newly compiled data set consists of 83 chloroplast genes sampled across 30 taxa that include chlorophytes and charophytes, including four members of the Zygnematophyceae, and land plants, that were sampled following a balanced representation of the main bryophyte and tracheophyte lineages. Analyses of non-synonymous site nucleotide data and amino acid translation data result in congruent phylogenetic trees showing the monophyly of bryophytes, with the Zygnematophyceae as the charophyte group most closely related to land plants. Analyses showing that bryophytes and tracheophytes evolved separately from a common terrestrial ancestor have profound implications for the way we understand the evolution of plant life cycles on land and how we interpret the early land plant fossil record.

The mitochondrial phylogeny of land plants shows support for Setaphyta under composition-heterogeneous substitution models.

Congruence among analyses of plant genomic data partitions (nuclear, chloroplast and mitochondrial) is a strong indicator of accuracy in plant molecular phylogenetics. Recent analyses of both nuclear and chloroplast genome data of land plants (embryophytes) have, controversially, been shown to support monophyly of both bryophytes (mosses, liverworts, and hornworts) and tracheophytes (lycopods, ferns, and seed plants), with mosses and liverworts forming the clade Setaphyta. However, relationships inferred from mitochondria are incongruent with these results, and typically indicate paraphyly of bryophytes with liverworts alone resolved as the earliest-branching land plant group. Here, we reconstruct the mitochondrial land plant phylogeny from a newly compiled data set. When among-lineage composition heterogeneity is accounted for in analyses of codon-degenerate nucleotide and amino acid data, the clade Setaphyta is recovered with high support, and hornworts are supported as the earliest-branching lineage of land plants. These new mitochondrial analyses demonstrate partial congruence with current hypotheses based on nuclear and chloroplast genome data, and provide further incentive for revision of how plants arose on land.

Phylogenomics provides robust support for a two-domains tree of life.

Hypotheses about the origin of eukaryotic cells are classically framed within the context of a universal 'tree of life' based on conserved core genes. Vigorous ongoing debate about eukaryote origins is based on assertions that the topology of the tree of life depends on the taxa included and the choice and quality of genomic data analysed. Here we have reanalysed the evidence underpinning those claims and apply more data to the question by using supertree and coalescent methods to interrogate >3,000 gene families in archaea and eukaryotes. We find that eukaryotes consistently originate from within the archaea in a two-domains tree when due consideration is given to the fit between model and data. Our analyses support a close relationship between eukaryotes and Asgard archaea and identify the Heimdallarchaeota as the current best candidate for the closest archaeal relatives of the eukaryotic nuclear lineage.

A meta-taxonomic investigation of the prokaryotic diversity of water bodies impacted by acid mine drainage from the São Domingos mine in southern Portugal.

The prokaryotic communities of water bodies contaminated by acid mine drainage from the São Domingos mining area in southern Portugal were analyzed using a meta-taxonomics approach with 16S rRNA gene sequences. Samples were collected in two seasonal sampling campaigns (summer and winter of 2017) from the most contaminated sites from where the water flows downstream to the freshwater reservoir of the river Chança. The physicochemical data indicate a trend of decreasing acid mine drainage contamination downstream of the mining area to the Chança's reservoir. The most contaminated sites (pH = 2.3-3.1) are distinguished by prokaryotic diversity with high abundances of operational taxonomics units related to acidophiles (genera Metallibacterium, Acidibacter, Leptospirillum, Acidobacterium, Thiomonas, Acidicapsa, Acidocella, Acidiphilium; family Acidobacteriaceae, order CPla-3 termite group). Likewise, in the transition zone in the mouth of the contaminated water flow into the Chança´s reservoir (pH = 6.4), a specific prokaryotic flora exists with some acidophiles, but notably with a cyanobacteria bloom and a high abundance of the genus Sediminibacterium (family I; order Subsection III). Moreover, the strong correlation between the abundance of acidophiles and characteristic physiochemical parameters (metals, acidity, and sulfate) confirm their potential as biomarkers of acid mine drainage pollution.

A haplotype-resolved draft genome of the European sardine (Sardina pilchardus).

BACKGROUND: The European sardine (Sardina pilchardus Walbaum, 1792) is culturally and economically important throughout its distribution. Monitoring studies of sardine populations report an alarming decrease in stocks due to overfishing and environmental change, which has resulted in historically low captures along the Iberian Atlantic coast. Important biological and ecological features such as population diversity, structure, and migratory patterns can be addressed with the development and use of genomics resources. FINDINGS: The genome of a single female individual was sequenced using Illumina HiSeq X Ten 10x Genomics linked reads, generating 113.8 gigabase pairs of data. Three draft genomes were assembled: 2 haploid genomes with a total size of 935 megabase pairs (N50 103 kilobase pairs) each, and a consensus genome of total size 950 megabase pairs (N50 97 kilobase pairs). The genome completeness assessment captured 84% of Actinopterygii Benchmarking Universal Single-Copy Orthologs. To obtain a more complete analysis, the transcriptomes of 11 tissues were sequenced to aid the functional annotation of the genome, resulting in 40,777 genes predicted. Variant calling on nearly half of the haplotype genome resulted in the identification of >2.3 million phased single-nucleotide polymorphisms with heterozygous loci. CONCLUSIONS: A draft genome was obtained, despite a high level of sequence repeats and heterozygosity, which are expected genome characteristics of a wild sardine. The reference sardine genome and respective variant data will be a cornerstone resource of ongoing population genomics studies to be integrated into future sardine stock assessment modelling to better manage this valuable resource.

Resolution of the ordinal phylogeny of mosses using targeted exons from organellar and nuclear genomes.

Mosses are a highly diverse lineage of land plants, whose diversification, spanning at least 400 million years, remains phylogenetically ambiguous due to the lack of fossils, massive early extinctions, late radiations, limited morphological variation, and conflicting signal among previously used markers. Here, we present phylogenetic reconstructions based on complete organellar exomes and a comparable set of nuclear genes for this major lineage of land plants. Our analysis of 142 species representing 29 of the 30 moss orders reveals that relative average rates of non-synonymous substitutions in nuclear versus plastid genes are much higher in mosses than in seed plants, consistent with the emerging concept of evolutionary dynamism in mosses. Our results highlight the evolutionary significance of taxa with reduced morphologies, shed light on the relative tempo and mechanisms underlying major cladogenic events, and suggest hypotheses for the relationships and delineation of moss orders.

Genomic blueprints of sponge-prokaryote symbiosis are shared by low abundant and cultivatable Alphaproteobacteria.

Marine sponges are early-branching, filter-feeding metazoans that usually host complex microbiomes comprised of several, currently uncultivatable symbiotic lineages. Here, we use a low-carbon based strategy to cultivate low-abundance bacteria from Spongia officinalis. This approach favoured the growth of Alphaproteobacteria strains in the genera Anderseniella, Erythrobacter, Labrenzia, Loktanella, Ruegeria, Sphingorhabdus, Tateyamaria and Pseudovibrio, besides two likely new genera in the Rhodobacteraceae family. Mapping of complete genomes against the metagenomes of S. officinalis, seawater, and sediments confirmed the rare status of all the above-mentioned lineages in the marine realm. Remarkably, this community of low-abundance Alphaproteobacteria possesses several genomic attributes common to dominant, presently uncultivatable sponge symbionts, potentially contributing to host fitness through detoxification mechanisms (e.g. heavy metal and metabolic waste removal, degradation of aromatic compounds), provision of essential vitamins (e.g. B6 and B12 biosynthesis), nutritional exchange (especially regarding the processing of organic sulphur and nitrogen) and chemical defence (e.g. polyketide and terpenoid biosynthesis). None of the studied taxa displayed signs of genome reduction, indicative of obligate mutualism. Instead, versatile nutrient metabolisms along with motility, chemotaxis, and tight-adherence capacities - also known to confer environmental hardiness - were inferred, underlying dual host-associated and free-living life strategies adopted by these diverse sponge-associated Alphaproteobacteria.

Improved phylogeny of brown algae Cystoseira (Fucales) from the Atlantic-Mediterranean region based on mitochondrial sequences.

Cystoseira is a common brown algal genus widely distributed throughout the Atlantic and Mediterranean regions whose taxonomical assignment of specimens is often hampered by intra- and interspecific morphological variability. In this study, three mitochondrial regions, namely cytochrome oxidase subunit 1 (COI), 23S rDNA (23S), and 23S-tRNAVal intergenic spacer (mt-spacer) were used to analyse the phylogenetic relationships of 22 Cystoseira taxa (n = 93 samples). A total of 135 sequences (48 from COI, 43 from 23S and 44 from mt-spacer) were newly generated and analysed together with Cystoseira sequences (9 COI, 31 23S and 35 mt-spacer) from other authors. Phylogenetic analysis of these three markers identified 3 well-resolved clades and also corroborated the polyphyletic nature of the genus. The resolution of Cystoseira taxa within the three clades improves significantly when the inclusion of specimens of related genera was minimized. COI and mt-spacer markers resolved the phylogeny of some of the Cystoseira taxa, such as the C. baccata, C. foeniculacea and C. usneoides. Furthermore, trends between phylogeny, embryonic development and available chemotaxonomic classifications were identified, showing that phylogenetic, chemical and morphological data should be taken into account to study the evolutionary relationships among the algae currently classified as Cystoseira. The resolution of Cystoseira macroalgae into three well supported clades achieved here is relevant for a more accurate isolation and identification of natural compounds and the implementation of conservation measures for target species.

Circulating small non-coding RNAs provide new insights into vitamin K nutrition and reproductive physiology in teleost fish.

BACKGROUND: Vitamin K (VK) is a fat-soluble vitamin known for its essential role in blood coagulation, but also on other biological processes (e.g. reproduction, brain and bone development) have been recently suggested. Nevertheless, the molecular mechanisms behind its particular function on reproduction are not yet fully understood. METHODS: The potential role of VK on reproduction through nutritional supplementation in Senegalese sole (Solea senegalensis) was assessed by gonadal maturation and 11-ketosterone, testosterone and estriol plasma levels when fed with control or VK supplemented (1250 mg kg-1 of VK1) diets along a six month trial. At the end, sperm production and quality (viability and DNA fragmentation) were evaluated. Circulating small non-coding RNAs (sncRNAs) in blood plasma from males were also studied through RNA-Seq. RESULTS: Fish fed with dietary VK supplementation had increased testosterone levels and lower sperm DNA fragmentation. SncRNAs from blood plasma were found differentially expressed when nutritional and sperm quality conditions were compared. PiR-675//676//4794//5462 and piR-74614 were found up-regulated in males fed with dietary VK supplementation. Let-7g, let-7e(18nt), let-7a-1, let-7a-3//7a-2//7a-1, let-7e(23nt) and piR-675//676//4794//5462 were found to be up-regulated and miR-146a and miR-146a-1//146a-2//146a-3 down-regulated when fish with low and high sperm DNA fragmentation were compared. Bioinformatic analyses of predicted mRNAs targeted by sncRNAs revealed the potential underlying pathways. CONCLUSIONS: VK supplementation improves fish gonad maturation and sperm quality, suggesting an unexpected and complex regulation of the nutritional status and reproductive performance through circulating sncRNAs. GENERAL SIGNIFICANCE: The use of circulating sncRNAs as reliable and less-invasive physiological biomarkers in fish nutrition and reproduction has been unveiled.